Decoherence induced by zero point fluctuations in quantum Brownian motion 
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We show a completely analytical approach to the decoherence induced by a zero temperature 
environment on a Brownian test particle. We consider an Omhic environment bilinearly coupled 
^ ' to an oscillator and compute the master equation. From diffusive coefficients, we evaluate the 

decoherence time for the usual quantum Brownian motion and also for an upside-down oscillator, 
, as a toy model of a quantum phase transition. 



(N 



PACS numbers: 03.65.Bz;03.70+k;05.40+j 



(N 

^ ■ I. INTRODUCTION 

in ■ 

I The emergence of classical behavior from a quantum system is a problem of interest in many branches 

■ of physics 0| . As it is well known, the quantum to classical transition involves two necessary and related 
conditions: correlations, i.e. the Wigner function of a quantum system should have a peak at the classical 

, trajectories; and decoherence, that is, there should be no interference between classical trajectories. To 

' study quantitatively the emergence of classicality, it is essential to consider the interaction of the system 

with its environment, since both, the loss of quantum coherence and the onset of classical correlations, 
^ • depend strongly on this interaction • Using this point of view, classicality is an emergent property of an 

\ open quantum system. The strength of the coupling between system and environment sets the decoherence 

. time which, roughly speaking, indicates the timescale after which the system can be considered classical 0]. 

I The very notion of quantum open system implies the appearance of dissipation and decoherence as an 

^ , . ubiquitous phenomena and plays important roles in different branches of physics (from quantum field theory, 

many body and molecular physics to theory of quantum information), biology and chemistry. Oftentimes, 
a large system can be described adequately as a composite system, consisting of two or a few subsystems 
(degrees of freedom) interacting with their environment (thermal bath) comprising a large number of degrees 
^ . , of freedom. Examples include electron transfer in solution, large biological molecules, vibrational relaxation 

■ of molecules in solution, excitons in semiconductors coupled to acoustic or optical phonon modes. Quantum 
processes in condensed phases are usually studied by focusing on a small subset of degrees of freedom and 
treating the rest as a bath. 

Decoherence is the main ingredient in order to find classicality. The interaction between the system and 
the environment induces a preferred basis which is stable against this interaction, and becomes a classical 
basis in the Hilbert space of the coupled system. Preferred pointer states are resilient to the entangling 
interaction with the bath. This "einselection" (environment induced superselection) of the preferred set of 
resilient pointer states is the essence of the environment. It is accepted that a rapid loss of coherence caused 
by the coupling with the environment is at the root of the non-observation of quantum superpositions of 
macroscopically different quantun states ■ A relevant property of the pointer states is their insensitivity 
to being monitored by the interaction with environment (and, therefore, are resistant to the entanglement 
caused by the environment). The less the states entangle, the more stable they are. All other states evolve 
into joint system-environment states, preserving their purity. 
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Our concern in this Letter is to analyze the effect of the zero point fluctuations of the environment, as 
a source of decoherence. The coupling of a quantum system to an environment generally leads to energy 
fluctuations in the test particle even at zero temperature Q . Since phases are time integrals of energy, zero 
point energy fluctuations make possible that decoherence occurs even at zero temperature. These fluctuations 
are a consequence of the finite coupling energy between the test particle (system) and the bath, and of the 
fact that the Halmiltonian of the isolated system does not commute with the interaction Hamiltonian. 

Vacuum fluctuations have several observable effects. The Lamb shift is a widely known example. Another 
one is the Casimir effect. In these examples the effect of vacuum can be thought in terms of the renor- 
malization of the original parameters characterizing the system. In contrast, the fluctuations we deal with 
in this Letter are not only absorbed into renormalized parameters of the test particle. Not only does the 
environment renormalize, but also it is a source of dissipation and noise for the system. Therefore, wc are 
considering the effect of quantum fluctuations of the environment over a quantum system, as the only source 
of decoherence (meaning that there is no possible thermal fluctuations inducing classicality) . 

The question about the influence of zero temperature environment on the interference phenomena has 
been discussed in the last years iTj . There have been studies on the temperature dependent weak local- 
ization measurements H , reporting residual decoherence in metals at zero temperature, in contradiction to 
theoretical predictions [9j , and on the zero-point decoherence induced by Coulomb interactions in disordered 
electron systems; just to mention a few examples. 

In previous works (see for example for an excellent review of the state of the art of decoherence) about 
decoherence in quantum Brownian motion, most of the conclusions are simply numerical or analytical only 
in the high temperature limit ^3 . Low temperature case was discussed in 0, 0| showing a numerical 
estimation of the decoherence rate. S. Sinha, in Ref. 0, studied the zero temperature case analytically. 
Under some approximations the author found an expression for the time dependence of the off-diagonal terms 
of the density matrix. In this Letter, we complete that study showing an exact calculation of diffusive terms, 
and also providing the decoherence timescale for different situations of interest. In addition, we solve the 
master equation for an upside-down Browni an p article to emphasize the role of zero temperature fluctuations 
during a second order phase transition |l3l Il4j . 



II. THE MASTER EQUATION AT T=0 



Let us consider a quantum particle (characterized by its mass M and its bare frequency fl) linearly coupled 
to an environment composed of an infinite set of harmonic oscillators (of mass to„ and frequency cij„). We 
may write the total action corresponding to the system-environment model as (we set h — I) 



S[x, q„ 



S[x] + S[q„] + Sint[x,q„ 



J^ds ^M{x^ -n^x^) + Y,lrnn{ql-ujlql) 



xq„ 



(1) 



where x and qn are the coordinates of the particle and the oscillators, respectively. The particle is coupled 
linearly to each oscillator with strength C„. 

The relevant objects to analyze the quantum to classical transition in this model are the reduced density 
matrix (obtained from the full density matrix integrating out all the degrees of freedom of the environment 
noted as g), and the associated Wigner function 



Pi-{x,x',t) 
W,{x,p,t) 



1 
2^ 



dq p{x,q,x\q,t) 

+ 00 



dye^Py p,(x + |,x-|,t). 



(2) 
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The reduced density matrix satisfies a master equation. Hu-Paz-Zhang [l5| have evaluated the master 
equation for the quantum Brownian motion problem (alternatively, one can write an equation of the Fokker- 
Planck type for the reduced Wigncr function in order to study the dynamics in phase space) 



p, = -i[H,y,t + ^MD?x\p,] -i-f{t)[x,{p,p,)] -D{t)[x, [x,p,]] - f{t)[x,[p,p,]]. 
The time dependent coefficients (in the case of weak coupling to the bath) are given by 



(3) 



lit) 
D{t) 

fit) 



M 



Mfl 



ds cos{D,s)t]{s) 



ds sin(ns)77(s) 



ds cos(f2s);/(s) 



(4) 
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ds sm(fls)iy{s), 



where 6Q^{t) is the shift in frequency, which produces the renormalizcd frequency that appears in the 
master equation, "/(t) is the dissipation coefficient related to the friction kernel defined below, and D{t) and 
f{t) are the diffusion coefficients, which produce the decoherence effects. Diffusion coefficients come from 
the noise kernel, source of stochastic forces in the associated Langevin equation. f{t) is named anomalous in 
the literature since it generates a second derivative term in the phase space representation of the evolution 
equation, just like the ordinary diffusion term f?(i) and 1^(1) wee the dissipation and noise kernels, 
respectively, 

poo 

Tjit) = / dLuI{Lu) sin ujt 
Jo 



doji [oj) coth — cos cut, 
2 



I{oj) is the spectral density of the environment, defined as I{uj) = (2/7r)A/7oA^ct;/((jj^ + A^) (where A is 
the physical high-frequency cutoff, which represents the highest frequency present in the environmct). In 
the high temperature limit of an ohmic environment (where litu) oc uj) the coefficients in Eq.Q) become 
constants. In particular, the diffusion coefficient can be approximated by £> ~ 2^QkBT M , where 70 is the 



dissipation coefficient [15j. In this limit, while 70 is a constant and D (x T, the coefficient / oc T~ can 
be neglected. Therefore, the term proportional to D is the relevant one in the master equation at high 
temperatures in order to evaluate decoherence. 

We will evaluate the time dependent coefficients of the master equation at zero temperature. For this, we 
set coth.phuj/2 = 1. These coefficients have been evaluated previously in Refs. [T^ll^ . 

The shift in frequency Sn^{t) is. 



performing integrations, we obtain 



duj / ds 
Jo 



a;2 + A2 



sin Ljs cos rJs, 



(5) 



sn'^it) -270 ^ 



A2 +f72 



1 - e 



-At 



cos fit r sin fit 

A 
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for times such that At > 1 the shift reads (see Fig. 1 (a)) 



= (7) 



A3 



'A2+172 

Dissipation coefficient (Fig. 1 (b)) comes ft-om the integral 



7(t) = A^ / dio I ds— — r- sinoissinJls, (8) 

TTii ./n ./o + A^ 



and it is given by 



l{t) = 70- 



A2 



'A2 + f72 

which has the foUowing asymptotic behavior 



1 



-At 



cos fit 



sin fit 



(9) 



7(0 = 70^^. (10) 
The normal diffusive coefficient (normally connected with decoherence effects) is coming from the integral 



Dit) = A / duj I ds— -T7 coswscosfls. (11) 

^ in in ^2 + A2 



This integral can be exactly solved. The result is: 



D{t) 



2M7n A2rj 

TT 172 + A2 



A 



Shi{At) ( — cos fit cosh At + sin fit sinh At 



Chi{Kt) cos nt sinh At + sin fit cosh At 
Si{m)] , 



(12) 



where, Chi{x) and Shi{x) are the hyperbolic Coslntegral and Sinlntegral respectively; Siix) is the Sinlnte- 
gral. 

The expression can be very well approximated, when At > 1, by 



2M70 A2f} 

This coefficient is the normal diffusion at T = 0. This is for any value of O. We are interested in time scales 
longer than the memory time 1/A. This coefficient is an oscillatory function of time. In fact, only in the limit 
fit ^ 1, as 5*1 goes to 7r/2, we can have an asymptotic value ~ M^qK^Q,/ [K"^ + fi^), independent of time 
|inl | (Fig. 1 (e)). In any other case, the coefficient has an initial transient and approaches the asymptotic 
value Doo as the 5*1 (see Fig. 1 (c)). It is important to note, that in the opposite case, when fit <C 1, 
normal diffusion is a linearly growing function of time, D ^ 2M {-^q / TT)K^n'^t / [A"^ + Vl?). similar to the result 
obtained in Ref.0- 
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FIG. 1: Time dependence of the coefficients of the master equation at T = 0. On top we show frequency renor- 
malization (a) and dissipation coefficients (b). Plots below show the normal D (c) and anomalous / (d) diffusion 
coefficients for short times in order to show the initial transient. Asymptotic values of diffusion are shown in (e) and 
(f). The parameters used in the plot are 70 = 0.05, A = fOO, Q. = 1. 



The anomalous diffusion coefficient is given by the integral 



7rS2 Jo 



It reads as 



f{t) = 270 



f72 + A2 

Chi{At) 



du) / ds— — ■ 



A 



cos LOS sin ris. 



Shi{At) ( ^ sin rJt cosh At - cos fit sinh At 



sin fit sinh At + cos fit cosh At 



Ci{nt) - log 



Again, for At > 1, this coefficient can be written as (Fig. 1 (d)) 



(14) 



(15) 



m = 2,o^^ (-c^m -io4 



(16) 



coefficient f{t) also approaches an asymptotic value when fit ^ 1, foo ~ — 27o(A2/(A2 + fl'^)) log A/ fl (Fig. 
1 (f)); and it docs to f{t) ^ — 27o(log At + F), when fit ^ 1 (F is the EulerGamma number). 

In Fig. 1, we show the time behavior of these coefficients. It is easy to probe asymptotic behavior also 
from analytical expressions. 
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It has been noted that the master equation in the high temperature hmit (or even in the 7 ^ constant 
approximation) has the pathology that the density matrix loses its positivity at short times (shorter than 
l/ksT). This violation is essentially due to the action of the friction term. Master equation Q does not 
have the pathological behaviour of the master equation at high temperature |15| . The dissipation coefficient 
is a time-dependent function that vanishes initially together with its first derivative. Therefore, the initial 
behaviour of the density matrix is diffusion dominated and positivity is preserved, even in the perturbative 
case, up to second order with respect to the coupling constant between system and environment. 

With these coefficients at hand, we will evaluate the decoherence time following Refs. [TollTl| . 



III. DECOHERENCE TIME AT T = 



We will analyze the decoherence process in a simple case. We prepare an initial superposition of delocalized 
(in position or momentum) states. We consider two wave packets symmetrically located in phase space, of 
the form [13: ^{x,t = 0) = ■^i{x) + ^2i x), where 

^1,2 = iVexp f- ^^2"^' ) cxp{±iPox), (17) 



l + exp(-^-<52p2 



(18) 



where is normalization, and 5 is the initial width of the wave packet. In terms of the Wigner function, 
the state at time t is W{x,p,t) = Wi{x,p,t) + W2{x,p,t) + Wint{x,p,t), where 



Si 



■ exp 



51 



exp {~5l{p TPc~ P{x =F Xc)f 



(19) 



and 



2N^ Si 



int 



(20) 



All the coefficients are functions of time, determined by the evolution propagator of the reduced density 



-S^ip - Pxf cos{2kpp + 2{kx - f3kp)x). 

n pre 

matrix and the initial state. The explicit form can be found in Ref. |lll |. The initial state is such that 
Sf = S2 = S'^, kx ~ Pq = Pci kp = Lq ^ Xc- kp and k^ indicate the evolution of the fringes in the momentum 
and coordinates directions of the phase space. 

As it was defined in the previous literature (see for example [Tol | ) , the effect of decoherence is produced 
by an exponential factor cxp(— ^int), defined as 



eM-An.) = I WU-^P)U^ ^. (21) 

^ [I^l(a;,p)|poakW2(x,p)| 

pcakj 

Initially, ^int = 0, and it is always bounded Ai^t < Lq/S^ + S^Pq = Ajntlmax- The fringe visibility factor 
Aint evolves in time as Amt = 4^D{t)kp — 4:f(t)kp{kx — (3kp). In the high temperature approximation, the 
anomalous coefficient is neglected and we obtain the very well known decoherence rate considering only the 
constant diffusion term, proportional to T . In our present case, at zero temperature, both coefficients D and 
/ contribute to the fringe visibility factor. A conservative choice is to assume fringes always stay more or 
less frozen at the initial values, and we can set kp = Lq and k^ ~ \/{2Lq). Neglecting the initial transient 
(i.e. Kt > 1), we use Eas. (|13|l and l|16(l to evaluate Ant- In order to have the simplest analytical expression 
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for the decoherence rate, we use a short time approximation to evaluate p, giving /3 ~ (see [ill). Thus, we 
get 

Ai^t-^LlD{t)~2f{t). (22) 

In order to evaluate the decoherence time tu, we have to solve 1 sa Aint{t — to)- From Ea. (|22|l it is not 
possible to find a global decoherence time-scale at T = 0. Nevertheless, we can find limits in which we are 
able to give different scales for decoherence. 

For example, for large natural frequency H, such as 51 ~ A {VLt ^ 1), it is easy to see that 

Aint ~ 2LlMjoAt + 470 (^t Ci{M) - , (23) 

giving a very short decoherence time-scale, 

\ — -. (24) 

This result will be valid as long as the product A/i§7o < 1, in order to be able to neglect the initial transient. 
It coincides with the decoherence time evaluated directly from Doc as in Ref.jlOf. In this limit, the anomalous 
coefficient docs not play any role (as we could check using /oo in H24(l ). 

In the opposite limit, when fit <C 1 (for times < t < ji)^ we can approximate Ajnt using asymptotic 
limit of 5*1 and Ci by 



8A^ 



ML. 



2 



27r 



^{ntf + t (logAi + r - 1) 



(25) 



resulting in a decoherence time bound, to < which could be large for very underdamped systems. Here, 
the logarithmic correction is due to the f{t) diffusion term (unlike Ref. 0, where anomalous diffusion was 
neglected). This scale is longer than the decoherence time in the high temperature, even in the case of 
low temperature for high natural frequency of the Brownian particle. Our result is still smaller than the 
saturation time tsat — 7(7^' time in which ^int reaches its maximum value. 

In the case we can neglect the second term in (|25|) [for example considering "macroscopic" trajectories 
(large MLq)], we can show 



Summarizing, in this Section we have shown analytical expressions for the decoherence rate at zero tem- 
perature. We were able to extract the decoherence timescales in different cases, giving new results respect 
to previous works 0, 0| and showing how to get known numerical results. 

IV. DECOHERENCE FOR THE UPSIDE-DOWN HARMONIC OSCILLATOR 

In this Section we are concerned with the analysis of the quantum to classical transition of the order 
parameter during a second order phase transition |l4l | . In a realistic model one should address this problem in 
the context of quantum field theory ^31 ■ This is a very difficult task since non Gaussian and non perturbative 
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effects are relevant. For this reason, we will only concentrate here in a toy model in ordinary quantum 
mechanics. 

Guth and Pi |l3l | considered an upside down harmonic oscillator as a toy model to describe the quantum 
behavior of this unstable system. This toy model should be a good approximation for the early time evolution 
of the phase transition, as long as one can neglect the non-linearities of the potential 0]. In this Section, 
we will analyze the decohcrence effects during a quantum phase transition in which the environment is at 
T = 0. 

Let us consider the unstable quantum particle (characterized by its mass M and its bare frequency f2) 
linearly coupled to a zero temperature environment composed of an infinite set of harmonic oscillators (of 
mass m„ and frequency ujn). As the coupling between system and environment is lineal, the result is exact, 
and can be easily obtained by replacing 17 by ifi in the Hu-Paz-Zhang equation. If the initial wave function 
is Gaussian, it will remain Gaussian for all times (with time dependent parameters that set its amplitude 
and spread). 

Let us solve Eqs.Q using a Gaussian ansatz for the reduced density matrix 

p,(E, A, t) = 7V(t)e-(2--C)A^g-(2a+C)S^g-4.bSA^ (27) 

while the reduced Wigner function is exactly evaluated as 

Wr{x,p,t) = -^r^e-(^'^+^)- e-^w, (28) 

where C{t) is a real function; and where S = l/2(x + x') and A = l/2(x — x'). 
The master equation, in the zero temperature limit, becomes 

a = 4:ab--/{t){2a~C) + D{t) + 2bf{t) 

1 

6 = -2(a2 - 2&2 - ^) - 2^^^ + ^ ^^"^W - (2a - C)f{t) 
C = ACb + 2{2a^C)-f{t)~2D{t)^4bf{t) 

N = 2Nb. (29) 
In this case, the temporal coefficients are given by (in the At > 1 limit) 



Ao2 27oA3 7oA2 



From Eas. l|27|l and H28|l we see that the relevant function to describe correlations and decohcrence is 
now 2a — C. For 2a — C = 0(1) wc have both correlations and decohcrence. The set of Eas. (|29|) can be 
solved numerically. In Fig. 3 wc show the behavior of 2a — C as a function of time. We see that it tends 
asymptotically to a constant of order one (of course the asymptotic value depends on the properties of the 
environment). 

The main conclusion of this Section is the following. In order to stud y a sudden quench quantum phase 
transition, at early times we can use the upside down potential 0,0,^3- When the system is isolated, 
due to the high squeezing of the initial wave packet, x and p become classically correlated Q. The density 
matrix is not diagonal. The "correlation time" depends on the shape of the potential. However, when 
the particle is coupled to an environment, a true quantum to classical transition takes place. The Wigner 
function becomes peaked around a classical trajectory and the density matrix diagonalizes. The decohcrence 



9 




0.2 0.4 0.6 O.e 1 

t 



FIG. 2; This coefficient measures the importance of the nondiagonal terms in the density matrix. Aiternativeiy, 
it is the width of the Wigner function. Is shows the rapid decolierence of the unstable particle coupled to a zero 
temperature environment. Parameters are 70 = 0.01, Q — 1, and A — 100. 

time at T = depends on the diffusion coefficients D and / and plays an important role in the early stages of 
a quantum phase transition, inducing classicality of the order parameter. Quantum aspect could be relevant 
if non-linearities are taken into account |l9j . Decohcrencc allows a classical description even in the nonlinear 
regime. 
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